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The rare decays B+ -> Df(p, B° -> D { * ] ~ K + and B + -> can occur only via 

annihilation type diagrams in the standard model. We calculate these decays in perturbative 
QCD approach. We found that the calculated branching ratio of B° — > K + agree with the 
data which had been observed in the KEK and SLAC B factories. The decay B + -> D ( * )+ K° 
has a very small branching ratio at O (10 -8 ), due to the suppression from CKM matrix 
elements. The branching ratio of B + — > Df(f> is of order 10 -7 which may be measured in 
the near future by KEK and SLAC B factories. The small branching ratios predicted in the 
standard model make these channel sensitive to new physics contributions. 



1 Introduction 

The generalized factorization approach has been applied to the theoretical treatment of non- 
leptonic B decays for years^. It is a great success in explaining many decay branching ratio!3. 
The factorization approach (FA) is a rather simple method. Some efforts have been made to 
improve their theoretical applicatioP and to understand the reason why the FA has gone weffl 
One of these methods is the perturbative QCD approach (PQCD), where we can calculate the 

annihilation diagrams as well as the factorizable and nonfactorizable diagrams. 

(*) 

The rare decays B -» DY'K(<j>) are pure annihilation type decays. In the usual FA, this 
decay picture is described as B meson annihilating into vacuum and the D { s * ] and K(<f>) mesons 
produced from vacuum then afterwards. To calculate these decays in the FA, one needs the 
Dg*^ — ► K((f>) form factor at very large time like momentum transfer O(Mb). However the 
form factor at such a large momentum transfer is not known in FA. The annihilation amplitude 



calculation of these decays is also unreliable. Here, we will try to use the PQCD approach, to 
evaluate the B — > D S K{4>) decays. By comparing the predictions with the experimental data, 
we can test the PQCD evaluation of the annihilation amplitude. 

A W boson exchange causes bd — ► cu or bu — ► d(s)c, which is usually described by the 
effective four quark operators, and the additional ss quarks included in D S K{4>) are produced 
from a gluon. This gluon attaches to any one of the quarks participating in the four quark 
operator. In the rest frame of B meson, s and s quarks included in D S K(4>) each has 0(Mg/4) 
momenta, and the gluon producing them has q 2 = 0(M^/4). This is a hard gluon. One can 
perturbatively treat the process where the four quark operator exchanges a hard gluon with ss 
quark pair. Therefore the quark picture becomes six-quark interactions. The decay amplitude 
is then expressed as product of the hard six quark operators and the non-perturbative meson 
wave functions. 



is a phenomenological parameter in QCD factorization approach 




QCDF 



2 Framework 



PQCD approach has been developed and applied in the non-leptonic B meson for 
some time. In this approach, the decay amplitude is separated into soft( ( I ) ), hard(i7), and 
harder(C) dynamics characterized by different scales. It is conceptually written as the convolu- 
tion, 

A~ Jd^hd^h Tr[C(t)$B(ki)$DSk2)$K(h)H(k 1 ,k 2 ,h,t)e~ s (% (1) 

where fcj's are momenta of light quarks included in each mesons, and Tr denotes the trace over 
Dirac and color indices. C(t) is Wilson coefficient of the four quark operator. In the above 
convolution, C(t) includes the harder dynamics at larger scale than Mb scale and describes the 
evolution of local 4- Fermi operators from mw, down to the scale t, where t = 0{\/KMb)- H 
describes the four quark operator and the spectator quark connected by a hard gluon whose q 2 
is at the order of t, and includes the hard dynamics characterized by the scale t. Therefore, 
this hard part H can be perturbatively calculated, which is process dependent. $m is the wave 
function which describes hadronization of the quark and anti-quark to the meson M. is 
independent of the specific processes. Determining <&m in some other decays, we can make 
quantitative predictions here. 

The large double logarithms (In 2 X{) on the longitudinal direction are summed by the thresh- 
old resummatioA'J, and they lead to St(xi) which smears the end-point singularities on X{. The 
last term, e~ s<yt \ contains two kinds of logarithms. One of the large logarithms is due to the 
renormalization of ultra-violet divergence In tb, the other is double logarithm In 2 b from the 
overlap of collinear and soft gluon corrections. This Sudakov form factor suppresses the soft 
dynamics effectively. Thus it makes perturbative calculation of the hard part H applicable at 
intermediate scale, i.e., Mb scale. 

In general, <&M,af3 having Dirac indices a, (3 are decomposed into 16 independent components, 

la/3) laBi (l5< J>1L ')ai3, (7^75)0/3, 75a/3- If the considered meson M is B or meson, to be 
pseudo-scalar and heavy meson, the structure (7^75)0/3 and 75^ components remain as leading 
contributions. Then, $m,o/3 is written by 

$M,at3 = -j= {{pMlh) a p<l>M + 75«/#m} • (2) 

As heavy quark effective theory leads to ~ Mb4>bi we have only one independent distribution 
amplitude for B meson. The heavy D s meson's wave function can also be derived similarly. 




Figure 1: Diagrams for B + — > Df<f> decay. The factorizable diagrams (a),(b), and non-factorizable (c), (d). 



In contrast to the B and D s mesons, for the K meson, being light meson, the 750"^ com- 
ponent remains. Then, K meson's wave function is parameterized as 



iSij 



75 Fz4>k{xzM) + m 0K j 5 (f)^(x 3 ,b 3 ) + m 0K >y 5 (/) jh 



(3) 



where mox = M\ / (m u +m s ) , v = (0, 1, t) , n = (1, 0, Ot) ■ InB — > D s (f> decay, only longitudinal 
polarization of the eft meson wave function is relevant, which is similar to K meson^. 

There are four kinds of Feynman diagrams contributing to the six quark hard dynamics, 
which is shown in Fig.l. t he ca lculation of the hard parts are tedious and channel dependent, 
the results are shown in refP^l 

3 Numerical evaluation 

In this section we show numerical results. First for the B meson's wave function, we use the 
same distribution amplitude as adopted in Ref.^l This choice of B meson's wave function is 
almost a best fit from the B — > Kit^tttt decays. For the Dl meson's wave function, we assume 
the same form as meson's on e? [ The wave functions <Px' P ' T of the K meson are expanded 
by Gegenbauer polynomials, which are given in Ref.^l 

For the neutral decay B° — > K + , the dominant contribution is the nonfactorizable 

annihilation diagrams, which is proportional to the Wilson coefficient ~ 1. The factoriz- 

able annihilation diagram contribution is proportional to a<i = C\ + C2/3, which is one order 

magnitude smaller. For the charged decay B + — > D^ + T^ '(</>) , it is the inverse situation. 

The propagators of inner quark and gluon in FIG. ^ are usually proportional to \jx{. One 
may suspect that these amplitudes are enhanced by the endpoint singularity around X{ ~ 0. 
However this is not the case in our calculation. First we introduce the transverse momentum 
of quark, such that the propagators become l/(xiXj + kj,). Secondly, the Sudakov form factor 
Exp[— S) suppresses the region of small k\. Therefore there is no singularity in our calculation. 
The dominant contribution is not from the endpoint of the wave function. As a proof, in our 
numerical calculations, for example, an expectation value of a s in the integration results in 
(a s /ir) = 0.10, Therefore, the perturbative calculations are self-consistent. 

The predicted branching ratios are^ 

Br(£° -> D~K + ) = (3.1 ± 1.0) x 10~ 5 , Br(B + -> D+K°) = (1.7 ± 0.4) x 10~ 8 , (4) 

Bt(B° -» D*-R + ) = (2.7 ± 0.6) x 10" 5 , Br(S+ -> D* s + K°) = (4.0 ± 0.8) x 10~ 8 , (5) 

for variation of the input para mete rs of wave functions. They agree with the experimental 
observation by Belli^ and BaBaP^, 

Bv(B° -> DJK + ) = (4.6l};? ± 1.3) x 10" 5 , Belle, (6) 
Br^ -> D~K + ) = (3.2 ± 1.0 ± 1.0) x 10~ 5 , BaBar, (7) 

and the experimental upper limit given at 90% confid enc e leveff^ Br(.B + — > Df7C°) < 1.1 x 
10" 3 . For B + -> D+4>, the predicted branching ratio is1^Br(5 + -> Dfcf)) = 3.0 x 10~ 7 , which 
is still far from the current experimental upper limit 15 : Bx(B+ -> D+(f>) < 3.2 x 10" 4 . 

Despite the calculated perturbative annihilation contributions, there is also hadronic picture 
for the B° — > D~ K + decay: B — ► D~ir + (p + ) — > D~ K + through final state interaction. Our 
numerical results show that the PQCD contribution to this decay is already enough to account 
for the experimental measurement. It implies that the soft final state interaction is not important 
in the B° — ► D~K + decay. This is consistent with the argument in Ref.^l We expect the same 
situation happens in other decay channels. 

4 Conclusion 



In two-body B decays, the final state mesons are moving very fast, since each of them carry more 
than 2 GeV energy. There is not enough time for them to exchange soft gluons. The soft final 



state intera ctio n may not be important. This is consistent with the argument based on color- 
transparencj 5 ^. The PQCD with Sudakov form factor is a self-consistent approach to describe 
the two-body B meson decays. Although the annihilation diagrams are suppressed comparing 
to other spectator diagrams, but their contributions are not negligible in PQCD approacH 

We calculate the B° D { * ] ~ K + and B + -> £>i* )+ if°(0) decays, which occur purely via 
annihilation type diagrams. T he bra nching ratio of B° — > D~K + decay is sizable, which has 
been observed in the B factoriePH^I The predicted branching ratio is in good agreement with 
the data. 
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